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Abstract 

Analytic expressions for the Fourier transforms of the Chebyshev 
and Legendre polynomials are derived, and the latter is used to find a 
new representation for the half-order Bessel functions. The numerical 
implementation of the so-called unified method in the interior of a convex 
polygon provides an example of the applicability of these analytic expres- 



1 Introduction 

The importance of Chebyshev and Legendre polynomials in numerical analysis 
is well known. The finite Fourier transform plays a crucial role in mathematics 
and applications. In spite of these facts it appears that the finite Fourier trans- 
form of Chebyshev and Legendre polynomials has not been constructed. Here, 
we present explicit formulae for these transforms. 

The finite Fourier transform of the Legendre polynomials can be expressed 
in terms of the half-order Bessel functions, thus a direct application of our 
results is the construction of an explicit representation for the half-order Bessel 
functions. 

Our results are motivated from the recent work of B Fornberg and collabo- 
rators [1|, |6| on the numerical implementation of the unified method of |2|-[4| 
to linear elliptic partial differential equations (PDEs) formulated in the interior 
of a convex polygon. In this case, the unified method yields a simple alge- 
braic equation, the so-called global relation, which couples the finite Fourier 
transforms of the given boundary data with the finite Fourier transforms of the 
unknown boundary values. For the determination of these boundary values 
one has to choose (a) appropriate basis functions and (b) suitable collocation 
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points in the Fourier space. Several such choices have appeared in the litera- 
ture [ 1 ], [5 1-[ 12 1; it appears that the best choice is (a) the unknown boundary 
values are expanded in terms of Legendre functions and (b) the collocation 
points are chosen on on rays introduced in 0, fl2l . We note that Chebyshev 
and Legendre functions give rise to equivalent finite bases of L 2 , and hence 
either choice will result in the same numerical method; however, in general 
the conditioning of the resulting linear systems will differ as conditioning is 
not invariant under matrix column operations. 

As an example of the applicability of these analytic formulae we apply this 
approach to the simplest possible polygon, namely to a square, and we use the 
explicit construction of the finite Fourier transform of the Legendre functions 
instead of the usual representation in terms of the Bessel functions. 

2 Fourier transforms of the Chebyshev and 
Legendre polynomials 

Theorem 1 Let T m (X) denote the finite Fourier transform of the Chebyshev 
polynomial T m (x), i.e. 



f m (X) = J ^ e- ax T m {x)dx, XeC, m = 0,1,2,..., 

where T m {x) denotes the Chebyshev polynomial 

T m {x) = cos(mcos _1 (x)), — 1<x<1, m = 0,l,2, 



Then, 



0, m= 1 



m 2 — 1 



, m = 0,2,3, 



Furthermore, 

m+l 

f m (A)=£c£ 

n=l 



-a 



(2.1) 



(2.2) 



(2.3) 



{ixy 1 " '< {ixy 

where the coefficients a'" are defined as follows: 

< = (-!)'", c^ = (-l) m+1 m 2 



AeC\{0}, m = 0,l,2, 



(2.4) 
(2.5) 



< = (-!) 



m—n+2 



r-^m £ 

k=l 



n + k-3 \"^r 3 



k-1 



j=k 

n = 3,4,..., m+l. (2.6) 
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Proof Making in (12. U the change of variables x = cos w, we find 



f m (X) = e- iUosw sinwcos(mw)dw, X e C, m = 0,l, 
Jo 

from which the results for X = follow. Using the identity 

•1 1 

— i A cosw 1 " 



sin we 



f g— iX cosvvj 

iX dw V / ' 

and employing integration by parts in (12.7b we find 
1 



(2.7) 



(2.8) 



iA / i \m 



f m (X) = tt e' A (-l 



where 



-mK„,(—iX) 



AeC\{0}, to = 0,1, 



(2.9) 



KJz)= e zcosw sm(mw)dw, z<E<C\{0}, m = 0,l,.... (2.10) 
Jo 

We now derive an analytic expression for K m (z). It is clear that Kq(z) = 0, and 
on replacing — iX by z in d2.81 > we find 



Ki(z)=- 



1 



= -(e z -e- z ), zeC\{0}. 



(2.11) 



The function K m (z) satisfies the relation 



K m+i (z) = K m {z)+K m _ x (z) + -{e? + {-\) m - l )e-:-), 

z z 



)m -i )e - z] 

zeC\{0}, m=l,2, 
For the derivation of this equation we will use the identity 

sin((m+ l)w) = sin((m— l)w) + 2 sin w cos (tow). 

Thus, 

Km+i (z) = r e zcosw sm{{m + \)w)dw 
Jo 

= -Km-i(z)+2 / sinwe zcos * cos(mw)dw. 



(2.12) 



(2.13) 



(2.14) 



Using integration by parts we find that the above integral can be rewritten in 
the following form: 



I sinwe ;:COSH 'cos(TOw)t/w 
Jo 



cos (tow) 



o Z Jo 



e ZCOSH 'sin(TOw)dw (2.15) 



1 -(e z + (-l)"^)e- z )--K m (z). 
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We now show that for m = 0, 1, . . ., 

K m {z) 

= ^ + (-1)^) 
1 m-1 /2\ n 

+ 7 I U) ((-l)V+(-l)V)£ 



n=l 



m— n / i r 1 

7 ) n 

(2.16) 



n+k-l 

n 



from which (12.41 ) follows using ( 12.9b . Clearly, this relation is valid for m = 
and m = I. 

Furthermore, the function K m (z) in ( 12.16b satisfies the relation ( 12.121 ) if, 
and only if, 



H=l 



1 m f1\ , "+7„±i._i\"+' ; - 1 



2m (2^ n 



E - ((-i)v+(-iroi 



-2 /o\ « 



-n-1 



ifc- 



m— n / ,7 1 



El7)((-i)V+(-i)"- 1 £ 



n+k-l 



"t-T'jnc-i-* 

(2.17) 



To prove ( 12.17b . we compare coefficients of each power of z on each side of 
this equation. The coefficient of z~ 2 on the left-hand side of ( 12.171 ) is given by 

m 

2{-e z + (- l) m+l e- z ) £ k(m + 1 - jfc) 



k=l 



_(^+(-iro 



i(m + 1 ) 2 m (m + 1 ) (2m + 1 ) 



(2.18) 



the coefficient of z 2 on the right-hand side of (12.171 ) is given by 



= -(^+(-l)» e -«, 



m—2 



2m 2 + 2££(m-l-&) 



2m 2 + (m- l) 2 (m-2)- -(m-2)(m- l)(2m-3) 



(2.19) 



and the expressions in ( 12.18b . ( 12. 19b are equivalent. The coefficient of z ( m+1 ) 
on the left-hand side of ( 12.17b is given by 



2"'((-l)'V + (-l)"' +1 e - z ) + l-;)=2 m ((-l)'V + (-l) m+1 e 



m+1 ^m!, 



;=i 



(2.20) 
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the coefficient of z ( m+1 ) on the right-hand side of ( 12.17b is given by 

m— 1 

- 2 m m{{- 1)"" V + {-l) m e-*) J] (« - ./) 

= 2 m m((-l)"V + (-l)" ,+1 e- z )(m- 1)!, (2.21) 

and the expressions in (12.201 i. ( 12.211 ) are equivalent. The coefficient of z~' n on 
the left-hand side of ( 12.171 ) is given by 



2 m -'((-l) m -V + (-l) m+, e- z )£ 



m + k — 2 N 



:2 m_1 (-l)" , - 1 (e r + «-*) 
:2 m m!(-l)" ,-1 (e z + Oi 



m— 1 



fj (m + 1 -./)+ m Y\ (m + l-j) 



(2.22) 



the coefficient of z m on the right-hand side of ( 12.171 ) is given by 



■2 m - 1 m((-l) m -V+(-l) m e^)£ 



2 / ,7 o \ m-\-k—3 
m + — 3 x 



t-7 ) n (»- 

./=* 



2 ffl - 1 m(-l) ffl -V + e - z ) 
2 m m(-l)'"- 1 (e z + e" z )(m-l)!, 



fe=l 

m— 2 m— 1 



rj(m-j)+('n-i)n( m -^ 



(2.23) 



and the expressions in ( 12.221 ). ( 12.231 ) are equivalent. Finally, we consider 
Z~( n+l \ n = 2, 3, . . . , m — 2. Omitting the common factor of 
2"((— l)"e z + (— l) m ~~ e~ z ) for convenience, the difference of the coefficients 
of the two sides of equation (12.17b is given by 



m—n+l 

E 

k=\ 



n+k-2 

n 



"t-f) no.-* 



n + k — 1 \ V-r , , A 



m—n— 1 

E 

/. = ! 



n + k-l 



n+k-l 



k _- ) n(— i-A 

(2.24) 



Furthermore, 



n + k-2 
k- 



n+k-2 



n + k-l 
k- 1 



n+k-l 



n+i-l 



./=/< 



1 n+k-2 



n + k-l 



- (m-k+1) 



J] (m-j). ( 2 - 25 ) 
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and 

m—n—l 

E 

k=l 



n + k-l 
k-l 



n+k-l 



n 



E 

<:=3 



fc-3 



11 (»-!-./) ( 2 - 2 6) 
"' " 1 7 B + jfc-1 \ (fe-l)(A:-2)(m-fe+l) n+fe - 2 



E 



t-i J ( B +jfc-i)( B +jfc-2) n ( m ^) 

Hence the above difference is given by 



m— n+1 

E 



n + ifc- 1 
fc-1 



n + fe- 1 



-(m-Jfc + 1) 1- 



(*-l)(*-2) 
{n + k- l)(n+k-2) 



n+Ar-2 



ri( w -i) 



m— ra+1 

« E 

^=1 



n + k-l 



k-l J n+k-l 



(n + 2k-3)(m-k+l) 
n + k-2 



n+k-2 



n ( m -j) 

j=k 



«E 

k=\ 



n + k \ k{n + 2k-m-l) n+ >^ 1 _ 



(2.27) 



where 



(m-k-iy.(n + k-2)Un + 2k-m-l) 

a k = y - LI // * =1, 2,..., m-n. (2.28) 

[m — n — k) ! [k — 1 J ! 

Hence a m -„+i-jt = — a*. If m — B is odd, then 

'm — n + 1 



n + 2 



m- 1=0, 



(2.29) 



and hence a,„-n+i = 0. Thus the coefficients of z ' n+1 ' of the two sides of 

— 2 — 

(12.17b are equal for n = 2, 3, . . . , m — 2. This establishes the validity of 02.171 1. 
which completes the proof of 02.4b . □ 
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Theorem 2 Let P m (X) denote the finite Fourier transform of the Legendre 
polynomial P m (x), i.e. 

p m {X) = J 1 e - iXx P m {x)dx, AeC, m = 0,l,2,..., (2.30) 

where P m (x) denotes the Legendre polynomial 
1 d m 

P m (x) = - - — ((x 2 -l) m ), -1< jc < 1, to = 0,1,2,.... (2.31) 



Then, 



Furthermore, 

m+l 

n=l 



^ 0l H 5; ::?; 2; .... (2.32) 



a -a 
_ — + (-l) m +"l — 
(ay y ' [iky 



, AeC\{0}, m = 0,1,2,..., 

(2.33) 



where the coefficients fi™ are defined as follows: 

J3j" = 1 //to even, (2.34a) 



7T7+7? — 3 \ / 777 +n — 3 



ft" = [(m + B)( + ^^(m), 



3,5,...,to+1 if m even, 
2,4,..., m+l if m odd, 

(2.34b) 



Rm / ^-1 \ ff / \ / 2,4,..., to ifmeven, 

n-l ) R ^+^-^ n= \l,3,...,mifmodd, 



(2.34c) 



where 



' 1, ,>t, V 7 1 0, 

(2.35) 

Proof We first consider ( 12.32b for to = 0: 

P (0) = f P (x)dx = 2. (2.36) 



Using the identity 



(2m + l)P m (x) = —(P m+1 (x)-P m -i(x)), -Kx<l, to = 1,2,..., 

(2.37) 

as well as the equations P m (\) = l and P m (— l) = (— l) m for m = 0,1,... we 
find the second line of ( 12.321 i: 

P m (0) = / P m (x)dx = ——[P m+X {x)-P m ^{x)] l _ x = 0. (2.38) 
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Form = 0, 1, P m (X) = T m (X) and thus P m {X) = f m (X); hence d2.41 > establishes 
the validity of d233l >. 

Furthermore, using (12.371 1 and integrating by parts we find 

4+i(A) = -^-(2m+l)P m (A)+4-i(A), AeC\{0}, m = l,2,.... 

(2.39) 

For convenience, we consider the following equivalent forms of (12.331 l for m 
even and odd respectively: 

m m — 1 

r~2j-'l )*mWu,,-;->(»') 



^) + ^-Z(A)«i, m -i(m), 



(2.40) 



and 

4(A) = -1 (i + zW - f( "-^ w+2) yW 

m-l '"r 3 

(-l)TiX (-1)' rt , n /m-;-l 



+ (r-2;- 1 i)^w)^+ 1 -^- 1 ( m )' 

(2.41) 

where 

Yr(A) = -{-(«-*-«?*) and = i^A)- (2.42) 
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For odd m, ( |2.39l > is equivalent to the following expression: 

m+l m ~ 1 



((2(m-i) + 3) ( m ™-/ + vr(A) + ( 



f(^)+ A J +1 ^I,m(m + 1)Z(A) 



(< 2 (»-^-<;-T- 2 0^)+(;-T- 2 ^ (A ») 



+ j?i,m-2(m-l)^(A) 
• p (2m + 1) ^1 + ^LlT ^i i( m )j #W 
~(2m+l)(ffl-l)(m + 2)^) 
V ; -(2m + l) 



A 2 

m — 3 
7=1 A 



((;ri" 1 i)*w+W"-/)+ i )(Viy 1 )' ,ti » 



(2.43) 



Furthermore, 



((-(-i)-i)(r_x_ 2 1 )^) 



m — j — 2 



m-2/-2 JZ(A) (m-1) 
-i ^ j ^ ' (2.44) 

Efc^(p(-^.)(;_T + i)"« 

m-j-l 

m - 2j 



)Rj, m -j-i(m-l). 
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It is straightforward to show that 

R j+ 1 ,m-j (m + l) = (2(m-j) + l)R j+ i, m -;-i ('«) , 

Rj,m-j- 1 (» - 1) = (2j + 1 )#,•+ 1 ,m-y- 1 (m) . (2.45) 

Hence for odd m, ( 12.391 ) is equivalent to the following equation: 

m+ 1 m ~ l 

^—j— I ^|(2(^i)+i)%,,,„- H w 

((2(-i) + 3) ( m * 2 7 +1 )vr(A)+ (™4.)z(A)) 

(-1)^ 

+ l A J +1 (2m+l) J R 1 , m _ 1 (m) Z (A) 

m+l 

= I ^|(2i + l)% 1 ,,»-MW 

3(-ir' 



1 m-1 



-«2,m-2(w)^(A) 



^2 ( 2m + 1 ) f 1 + ^lZ R hm- 1 (m)j Z (A ) 
^ ( 2m + 1) (m - 1) (m + 2)i/a(A) 
^ (2m+l) 

m-3 

;=i A 



((;rt- 1 i)^ i>+(2(m - J)+1, (T4 1 )^) 



(2.46) 
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Hence, 

(-1)^ 
A 

m—3 

' E T^27 /? ./+i.<»-./-iH 



(2(» - i) + 1) ((2(» - J) + 3)( m ™ "/ + J - (2; + l)^-; 

-(2» + l)(" B ^- 1 )) ¥ r W 

(-I)V 

+ A „,_! /g 2 ,m-2(w) [(2m - 1 ) ( (2m + 1 ) V (A ) + (m - 1 )j (A ) ) 

-(2m+l)((m-2)^(A) + (2m-l)vr(A))] 
- ^ ((m + 2) Q (m + 4) (m + 1 ) (m - 1 ) y(A) + 1 (m + 1 )x (A) 

-m^(m + 2)(m-l)(m-3)vr(A) + i(w- 1 l*(A) 

/(-l)^ 3(-l)^ 
+ 1 Am+ i (2w + l)*i,™-i(m) + Am li ^2,m-2(m) 

+ ^ (2m + 1 ) (l + { -j^rR 1 ,m- 1 (m)^ % ( A ) 

+^a(2m+l)(m-l)(m + 2)vr(A) 
= 0. 

Equation ( I2.47l i follows from the identities 

(2(-i) + 3)( m !'-/ +1 )-(2/ +1 )(;4.;> 1 ) 

-(2 m + i)( m m -j''- 1 ) 0. ,:.4S, 



(2(m-;) + l) -(2/ + 1 )' 



(2.47) 



m — 2j J \ m — 2 j 
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(-1)^ (-1)^ 

and the coefficients of ^ m+1 R\, m -\{m)x{X), t \l/(X)R 2 , m -2(m), 

^ ol-i" %{^) R l,m-z{m), ip(A) and -^pZW on the left hand side 



of (2 .471 > vanishing: 

(2m + 1) - (2m + 1) = 0, (2m - l)(2m+ 1) - (2m + l)(2m - 1) = 0, 
(2m - l)(m - 1) - (2m + l)(m - 2) - 3 = 0, 
(m - 1 ) (m + 2) ((m + 1 ) (m + 4) - m(m - 3) - 4(2m + 1 )) = 0, 
(m + 1) (m + 2) - m(m - 1) - 2(2m + 1) = 0. 

The validity of ( 12.391 1 for even m can be demonstrated similarly, completing 
the proof of equation ( 12.331 1. □ 

3 An explicit representation for the half-order 
Bessel functions 

Lemma 3 Let J m+ i (A) denote the half-order Bessel function, i.e. 
J i (A) = — / cos ( ( m H — ) T — AsinT)dT 

m+ 2 V ' %h \\ 2 J ) 

(_l)m+\ /•<*>/ ( l\ \ (3 1) 

+ - — ^ J expl -AsinhT- I m + - J r jdx, 

AeC, m = 0,l,2,.... 

Then 

4+'( ) = > m = 0,1,2,..., (3.2) 
and J m+ i (A) admits the following explicit representation: 



1 m+1 

j jm = -= r /3™ 



-iX 

j \ n+m 



AeC\{0}, m = 0,1,2,..., (3.3) 

where the coefficients j3™ are defined by equations (12.34b . 

Proof Recall that the finite Fourier transforms of the Legendre polynomials 
P m {k) can be expressed in the form 

UV = ±;]ffj m+k (V- (3-4) 

This equation, together with the results of Theorem[2] immediately gives equa- 
tion @3). □ 
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Figure 1: The square with sides of length 2. 



4 The modified Helmholtz equation in the inte- 
rior of a square 

Let u(x,y) satisfy the modified Helmholtz equation 

d 2 u d 2 u 



dx 2 dy 2 



4u = 0, 



(4.1) 



in the interior of a convex polygon. Associated with each side {Si}" of this 
n-gon, define the following functions 



Uj{X) - 

and 

Uj{X) - 



"My 



A + — I u 



dx- 



ik + j^- I u 



dy 



j=l,2,...,n, AeC\{0}, (4.2) 



■u, 1 A + — i a 



dx - 



ik + — i // 



dy 



7 = l,2,...,n, AeC\{0}. (4.3) 
Then, the functions {uj(X)}" and {"/(A)}" satisfy the global relations lfl2l 



£«;(A) = 0, £fiy(A) = 0, AeC\{0}. 



(4.4) 



For brevity of presentation, we consider the simplest possible polygon, namely 
the square with corners at (Figure Q]) 

(-1,1), (-1,1), (1,-1), (1,1). 
For the sides Si, S%, S3, S4, we have respectively 

z = — 1 + z = x — 2, z = l + £y, z = x + /. (4.5) 
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Hence, taking into account the orientations of the sides, we find the following 
expressions: 

-l 



J+l 



(4) 

"My - 



dy, (4.6a) 

<fo, (4.6b) 

dy, (4.6c) 

c/jc. (4.6d) 



Let Dj and denote the parts of uj corresponding to the Dirichlet and Neu- 
mann boundary values. Then, 

u l {X) = -e( iX+ ^)^ 1 (X)-(iX + ^-)e( a+ -k)£> 1 {X), (4.7a) 

a 2 (X) = -e-( l+ i)fi 2 (-iX) + (x + j^j e-( l+ j)D 2 (-iX), (4.7b) 
u 3 (X) = e-( a+ K)ft 3 (X)+ (iX + -^e-( iX+ rx)£> 3 (X), 
w 4 (A) = e( X+ x)N 4 {-iX) -U + jj e( X+ x)b 4 {-iX). 



(4.7c) 
(4.7d) 



For simplicity, we consider the symmetric Dirichlet boundary value problem: 



(— l,y) = cosh(l)cosh(v / 3y)+cosh(v / 3)cosh(3'), — 1 <y< 1; 



(4.8a) 



: u(l,y) = cosh(l)cosh(v / 3y)+cosh(v / 3)cosh(y), — 1 <y< 1; 



(4.8b) 

ift> = u(x, — 1)= cosh(l)cosh(v / 3x)+cosh(v / 3)cosh(x), — 1 <x< T, 

(4.8c) 

m^ = m(jc, 1) = cosh(l)cosh(v / 3-'c)+cosh(v / 3)cosh(x), — 1<jc<1. 

(4.8d) 



Then. 



Thus. 



u(x,y) =u(—x,y), u(x,y) — u(x,—y), u(x,y) = u(y,x). 



P) 



(1) (4) 

U X , Uy 



(2) (2) _ (1) 

Uy , Uy Ux 



(4.9) 
(4.10) 



Hence, the first of the global relations ( 14.41 i becomes 



cos U - - Ni(X) +cos /A - — Ni(-iX) 



iX 



sin A- T )£> 1 (X)+ [iX 



iX 



iX 



iX 



£>i(-iX), 



AeC\{0}. (4.11a) 
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The simplest way to obtain the second of the global relations ( 14.4b is to take 
the Schwartz conjugate of equation d4.11al > (i.e. to take the complex conjugate 
of (14.1 lab and then to replace X with X). This yields the equation 

cos (x - ij #i(A) + cos \iX - -^J #1 (iX) 

= (A - I) sin (A - 1) A W + (/A - 1) sin (/A - £) Di (/A), 

AeC\{0}. (4.11b) 

Using N basis functions to approximate u x l \ equations d4.1 U yield 2 equa- 
tions for N unknowns. 

Regarding the numerical solution of the global relations ( 14.1 U . Fourier 
basis functions Q, lT7l- lfT2l . as well as Chebyshev and Legendre polynomials 

|0D> ED, have been used to approximate u{ l \ In most of the earlier papers 
the collocation points A G C \ {0} were chosen to lie on the rays in the com- 
plex A -plane which are parallel to the edges of the polygon and its reflection 
in the imaginary axis. Recently B Fornberg and collaborators introduced the 
use of the so-called Halton nodes, [1 1, [6|. 

It appears that the most efficient numerical method involves (a) approxi- 
mating the data Ux ' in terms of Legendre polynomials (following Fornberg) 
and (b) using the collocation points employed in our earlier work lfl2l where 
ideas of A G Sifalakis and collaborators [9| for the Laplace equation were ex- 
tended to the modified Helmholtz equation. Furthermore, in order to ensure 
that the collocation matrix remains well conditioned as the number of basis 
functions increases, it is important following Fornberg to (i) divide each row, 
as well as each column, of the collocation matrix by its ^-norm |1| and (ii) 
to "over-determine" the linear system by choosing the number of collocation 
points to be about the same as the number of unknowns. 

Numerical experiments suggest that Legendre polynomials yield spectral 
accuracy rather than the algebraic accuracy found with a Fourier basis. Fur- 
thermore, as a result of choosing the collocation points to be on the above 
rays, the semi-block circulant structure of the collocation matrix for regular 
polygons, demonstrated for the Laplace equation in |8|, is preserved in modi- 
fied Helmholtz equation as well. 

Plots of the relative error (defined in [ 12|), as well as of the matrix con- 
dition number as a function of N, for N/2, N, 3N/2 and 2N collocation points, 
are presented in Figure [2] The rectangular collocation matrix was inverted by 
using the "backslash" command in Matlab. It is clear that over-determining 
the linear system by a factor of 2 is sufficient to achieve very good matrix 
conditioning. 
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Figure 2: Numerical solution of the global relations ( 14. 1U for the symmetric 
Dirichlet boundary value problem (14.81 ). 
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